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The second moment of the pion's distribution amphtude * 
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We present preliminary results for the second moment of the pion's distribution amplitude. The lattice for- 
mulation and the phenomenological implications are briefly reviewed, with special emphasis on some subtleties 
that arise when the Lorentz group is replaced by the hypercubic group. Having analysed more than half of the 
available conflgurations, the result obtained is — 0.06 ± 0.02. 



1. INTRODUCTION 

The distribution amplitude 0(a;i,a;2) gives the 
probability for finding two coUinear partons with 
fractions xi and X2 of the meson's momentum [Q. 
It plays an important role e.g. in exclusive hard 
scattering processes and in non-leptonic decays 
of heavy mesons, since it allows one to separate 
a hard scattering amplitude, computed in pertur- 
bation theory, from the effect of the soft Fourier 
modes, whose dynamics is non-perturbative. 

For example, the relevant diagram for the pion 
electromagnetic form factor is shown in Fig. ^ 
The form factor is defined as: 




Figure 1. Photon-pion interaction, factorised into 
a hard (perturbative) scattering amplitude and a 
non-perturbative part described by the pion dis- 
tribution amplitude. 
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and can be written in terms of the distribution 
amplitude: 



can be extracted from the matrix element of the 
lowest-twist operator appearing in the OPE of 
Eq.|: 



[dx][dy] ^\x, Q^)Th{x, y; Q^)cl>{y, Q^) (2) 



where [dx] = da;ida::2(5(a;i -|-a;2 — 1), and Th is the 
perturbative amplitude computed at the quark 
level. 

The amplitude is given by the non-local 
Fourier transform of the matrix element of the 
fermionic bilinear lEl: 



where 
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<l>tp{x, Q2) ^ FT [(0|rC(^i)^^(^2)k)] 
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A light-cone expansion relates (p to fermion bilin- 
ears with covariant derivatives. In particular the 
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and the ellipses indicate terms that are propor- 
tional to . 

The distribution amplitude has been studied in 
recent years using sum rules (see results quoted 
in 0]) and lattice simulations [§,1). The accu- 
racy of the lattice determinations so far has not 
been good enough to compare precisely with sum 
rules predictions. The aim of the current study is 
to use modern lattice technology to improve the 
precision of the result. 
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2. LATTICE DETAILS 

The following UKQCD quenched configura- 
tions have been used to compute the relevant ma- 
trix elements. 

• (3 = 6.2, a"i = 2.64 ± 0.10 GeV, 24^ x 48 
lattice; 

• SW action, with Csw — 1.61; 

• the three k values used in our simulation 
(0.13460, 0.13510, 0.13530) correspond to 
physical light pseudoscalar masses ranging 
from 350 to 850 McV. 

• 180 configurations are available, the results 
presented here are based on the analysis of 
a subset of 129; 

• non-improved local operators have been 
used, without fuzzing of the light quarks. 

3. POWER-LIKE DIVERGENCIES 

In this work, we concentrate on the second mo- 
ment of the distribution amplitude, i.e. to = 2 in 
Eq. H and ^ above. The choice of the Lorentz in- 
dices for 0/ioMiM2 is crucial in order to avoid mix- 
ing with lower-dimensional operators. We choose 
jj^o to be the time-direction. It is convenient, but 
not necessary, to avoid time-derivatives and hence 
we do not symmetrise /jq with /Xi,/i2- Using the 
notation of Q to classify the irreducible represen- 
tations of the hypercubic group, we obtain: 

• Op[^cr], with [i^v^a^[i and symmetrised 

over V and cr, transforms like a (i, i) © 8 
reducible representation, and therefore does 
not mix with lower dimensional operators. 
However the 8 irrep leads to a term propor- 



where r is the 
which needs to be 



tional to 

fourth available index 
subtracted in order to have an operator pro- 
portional to ppp^pcr, as in the continuum 
limit. This subtraction is performed in the 
definition of R\ below. Note that two spa- 
tial components of p need to be non-zero to 
get a non-vanishing signal. 



0^lvv, with [i^v transforms like a (^, ^)©8 
representation; however the subtracted op- 
erator: 
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simply transforms like 8. In this case there 
is no mixing, and a non-zero signal is ob- 
tained with just one non-vanishing compo- 
nent of p = (1,0,0). 

4. PRELIMINARY RESULTS 

The relevant matrix elements are extracted 
from suitably defined two-point functions. For 
a generic operator Q: 
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where Z = (7r(p)|^/'75^/'(0)|0). 
For large time separations: 
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Only the results for the heavier k are presented 
here. The pion propagator is shown in Fig. ^. 
The effective mass and single-exponential fits give 
consistent results for 6 < t < 15. Since the rel- 
ative error on the signal of interest for the sec- 
ond moment becomes large for t > 15, the fit- 
ting ranges chosen in this work are shorter than 
the ones usually used for light spectroscopy. De- 
spite this discrepancy, our value for /^r agrees with 
other UKQCD determinations from the same set 
of configurations within 2a. 

The time dependence of i?2, obtained from 129 
configurations, is given in Fig. |. We have cur- 
rently analysed Ri on 55 configurations only. In 
order to compare the two determinations, the re- 
sults for Ri are also displayed in the same figure. 
The agreement between the two sets of results is 
good. At this preliminary stage R2 alone is used 
to extract our final number. 
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Figure 2. Pion propagator from 129 configura- 
tions at K = 0.13460 
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Figure 3. Ratios of two-point functions i?i and 
R2 for K = 0.13460. The results for Ri and R2 
are obtained from 55 and 129 configurations re- 
spectively. 

5. CONCLUSIONS 

Fitting the points in Fig. || to a constant yields 
(^2^^ ^ 0.06 ±0.02 

for the lattice value of the second moment of the 
pion distribution amplitude. The analysis of the 
complete set of available configurations for both 
Ri and R2, together with a full discussion of the 
systematic errors, will be presented in a forth- 
coming publication The statistical error in 
our result will of course decrease as the num- 
ber of configurations is increased. The multi- 



plicative renormalisation of the operator has yet 
to be calculated, either perturbatively or non- 
pert urbatively. 

Our preliminary result for the value of {^'^) sug- 
gests that it is lower than the value predicted 
by sum rules. Such a result could be explained 
for example by a more peaked distribution for ^ 
around the origin. 

Computations of the distribution amplitudes of 
the p meson and the proton are currently un- 
der way. Our analysis can also be extended to 
the case of non-degenerate quarks and so a study 
of the distribution amplitudes of the K and B 
mesons is feasible. 
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